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Abstract 

Spectral analyses of static pressure fluctua- 
tions measured in turbine engine combustors at low 
engine speed show good agreement with theory. At 
idle speed the high pressure turbine is unchoked. 
Above idle speed the turbine chokes and a signifi- 
cant change in the shape of the measured combustor 
pressure spectrum is observed. A simplified theo- 
retical model of the acoustic pressure generated 
in the combustor due to the turbulence-flame front 
interaction did not account for acoustic waves 
reflected from the turbine. By retaining this 
simplified combustion noise source model and adding 
a partial reflecting plane at the turbine and com- 
bustor inlet, a simple theoretical model has been 
developed that reproduces the undulations in the 
combustor fluctuating pressure spectra. Plots of 
the theoretical combustor fluctuating pressure 
spectra are compared to the measured pressure 
spectra obtained from the CF6-50 turbofan engine 
over a range of engine operating speeds. Reason- 
able agreement exists, and it is thus concluded 
that the simplified combustion noise theory when 
modified by a simple turbine reflecting plane ade- 
quately accounts for the changes in measured com- 
bustor pressure spectra. It is further concluded 
that the shape of the pressure spectra downstream 
of the turbine, neglecting noise generated by the 
turbine itself, will be the combustion noise spec- 
tra essentially unchanged except for the level 
reduction due to the energy blocked by the turbine. 

Introduction 

The reduction of aircraft noise through the 
use of high-bypass-ratio turbofan engines has 
revealed another noise source for concern: core 

engine noise. Core engine noise is defined as the 
noise produced by the gas generator-turbine- 
tailpipe combination used to provide shaft power to 
drive the fan in a turbofan engine. Reference 1 
summarizes various combustion noise prediction 
methods for aircraft engines available before 1974. 
The methods include both engine and combustor rig 
predictions. 

More recently it has been recognized that gas 
generators in earth to orbit propulsion systems 
create fluctuating pressures that load structures 
and contribute to high cycle fatigue failure of 
critical structures. When prediction of the gas 
generator fluctuating pressure is required to 
determine high cycle loading the methods used to 
calculate core engine noise can be applied. 

Since the publication of Ref. 1 a data bank 
from various combustion noise experiments has been 
created by NASA. Additionally, a number of 
engine tests have been performed 10-19 to determine 
the amount of core engine noise reaching the far- 
field observer. The acoustic power determined from 
the data in Refs. 2 to 9 has been correlated with 


engine operating parameters in Ref. 20. A number 
of theoretical investigations is reported in the 
literature: Ref. 21 sets forth several possible 
combustion noise source terms to be evaluated; 

Ref. 22 gives a correlation of data, based on 
theory, for the acoustic power in a combustor with 
engine operating parameters; Ref. 23 gives scaling 
laws for combustion noise radiating from open 
flames that are based on a theory evolved by using 
the chemical reactions; Ref. 24 calculates the 
1/3-octave spectrum for a propane-air open flame; 
and Ref. 25 points out the need to know the heat 
release distribution in the combustor before the 
dynamic design can be performed. References 21 to 
24 have focused attention on combustion noise 
theory and have presented a fundamental, though at 
times complex, picture of possible combustion noise 
mechanisms. Reference 26 has shown that the com- 
bustor fluctuating temperatures and pressures are 
correlated in ducted burners, that heat release 
distribution is important and that the fluctuating 
volumetric heat release rate needs to be investi- 
gated as a source term in combustion qenerated 
noise. 

In an effort to understand the dominant com- 
bustion noise mechanism and to provide a theoreti- 
cally based prediction that allows insight into the 
noise-generating parameters in ducted combustors, 
a simple combustion noise model was derived in 
Ref. 27. The fluctuating pressure spectrum gener- 
ated by the interaction of the turbulence with the 
mean energy gradient in a burner connected to in- 
finite ducts (no reflected waves in either direc- 
tion) was determined analytically and compared to 
pressures measured in a CF6-50 turbofan engine com- 
bustor. While good agreement existed between data 
and the theory, deviations existed above engine 
idle condition that need to be explained. In the 
derivation of the acoustic source pressure in 
Ref. 27 the effect of acoustic waves reflected 
from both the turbine and combustor inlet planes 
was neglected. 

The objective of the present work is to pro- 
vide a theoretical solution for the acoustic pres- 
sure generated by a turbulent-flame interaction 
noise source located in a duct between two reflect- 
ing planes and to demonstrate that accounting for 
simple reflections can explain the variation in 
the narrowband spectral data from the existing 
theory of Ref. 27. All symbols used herein are 
defined in Appendix A. 

Theoretical Model for Combustion Noise with 
Ref lectinq Planes' 

The model for combustion noise considers a 
one-dimensional constant area, duct flow (Fig. 1). 
The flow is composed of a mixture of fuel and air. 
Turbulence is generated upstream of the flame front 
by turbulators. The flame front is stabilized 
aerodynamical ly at a point just downstream of the 
turbulators but does not touch the turbulators. 
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The source region of combustion-related noise 
starts at the flame front. To simplify the calcu- 
lations, the flame front is considered to be in a 
plane perpendicular to the flow direction. The 
regions upstream and downstream of the source 
region have constant though different properties 
and are terminated with partial reflecting planes. 

In engine combustors the combustion zone tempera- 
ture varies with axial distance. In this analysis 
the wave number is assumed constant in the region 
downstream of the flame front thus neglecting 
variation in sonic velocity in the combustion zone, 
he reflecting planes represent the turbine and com- 
bustor inlets in a real turbine engine or turbopump. 

Infinite Duct Theory 
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The narrowband fluctuating pressure spectrum 
measured in the combustor of the CF6-50 engine' has 
been predicted using a simple source theory in an 
infinite duct. 2 ' A comparison of the theory with 
data is shown in Fig. 2 for the engine operating at 
both low and high thrust. At the idle condition 
(Fig. 2(a)) the theory and data agree very well if 
the closely spaced minimum value points in the 
spectrum are ignored. At the 99.8 percent thrust 
condition (Fig. 2(b)) deviation of the measured 
spectra from that predicted from the infinite duct 
theory is apparent at frequencies around 300 and 
1300 Hz. To account for this deviation a theoreti- 
cal analysis including the effects of reflected 
waves on the measured spectra has been developed. 
The model, as expected, does not explain the 
closely spaced minimum value dips in the spectra. 

Reflecting Plane Theory 
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This equation was programmed and integrated numeri- 
cally over the combustion zone using a computer. 

The output of the program was the power spectral 
density of the acoustic pressure. 


The one-dimensional Fourier-transformed acous- 
tic wave equation with the turbulence-flame front 
noise source term from Ref. 27 is written here as: 


For comparison purposes the acoustic pressure 
calculated for the same combustion noise source but 
in an infinite duct is, from Ref. 27: 
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For low Mach number flow in the combustor the right 

side of equation (1) (source forcing function) 
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the solution of equation (1) with source term ex- 
pression (2) is accomplished using the Greens func- 
tion. First the duct is divided into two regions: 
one, upstream of a point source located in the 
source region and the other downstream of the point 
source. Complex reflection factors are defined and 
the method used for solving the wave equation for 
a string with source and reflecting points is fol- 
lowed as given in Ref. 28 (page 132-138) (see 
Appendix B for the derivation). The equation for 
the acoustic pressure in the downstream region is 
from Appendix B (Eqs. (B15) and ( B16 ) ) : 


(5) 

It can be shown that for equal wave numbers 
(ky = k 0 ) and for both reflection factors and phase 
angles equal, zero, equation (3) reduces exactly to. 
equation (5) as might be expected. 

Trends Predicted From Theory 

The acoustic pressure as a function of fre- 
quency has been calculated using Eqs. (3) and (4). 
This calculation accounts for reflections from the 
turbine and the combustor inlet. An acoustic pres- 
sure efficiency using infinite tube theory, was 
calculated in Ref. 27 using fluctuating pressure 
data measured in a CF6-50. turbofan combustor at the 
3.8 percent thrust level.' 
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The acoustic pressure efficiency as determined 
in Ref. 27 was 0.045. However, due to the omission 
of a radical in the equation used to calculate the 
acoustic pressure efficiency in Ref. 27 and because 
of terms in the frequency parameter, a corrected 
value of 0.088 has been used with the present theo- 
retical calculations. The narrowband data,' used 
herein were analyzed using a two Hertz bandwidth 
filter. All theoretical pressure spectra for com- 
parison purposes have been multiplied by the two 
Hertz bandwidth so that they are directly compar- 
able with the as analyzed data reported in Ref. 7. 

Measurement Location Effects 

The measured fluctuating pressure spectrum 
reported in Ref. 7 was obtained at the borescope 
opening in the combustor just downstream of the 
combustor inlet station. Figure 3 shows the theo- 
retical pressure spectrum at 100 percent thrust 
for duct measurement locations x/i = 0.1 through 
0.4 with the same reflection factor (R = 0.35) and 
phase shift (» = n/2) at both reflecting planes. 
Little difference exists below 1000 Hz over the 
range of axial locations shown (see Fig. 3); how- 
ever, significant difference exists above 10Q0 Hz. 

Reflection Factor Effect 

The magnitude of the reflection factor has a 
large effect on the spectral shape as shown in 
Fig. 4 for the CF6-50 combustor operating at 100 
percent thrust. Figure 4(a) presents the theoreti- 
cal narrowband pressure spectrum at x/i = 0.2 for 
equal reflection factors from R = 0.0 to 0.8 and 
a phase angle, * = n/2 radians. Large changes of 
the spectral shape are noted with changes of R in 
this range. For complete reflection, (Fig. 4(b)), 
Rj = Rg = 1, the resonant frequency calculated from 
f = C/4i for the upstream and downstream sonic 
velocities respectively are on the order of 386 and 
529 Hz where the values are calculated by treating 
the combustor as if it were entirely cold or en- 
tirely hot respectively. Within this range of fre- 
quencies several oscillations occur in the spectra. 
As will be shown later the lower values of the 
reflection factor appear to be more representative 
of the data. 

Variation of the phase angle associated with 
the reflection factor, as shown in Fig. 5, has a 
significant effect on the narrowband spectral 
shape. Figure 5 presents the theoretical pressure 
spectrum at the CF6-50 engine 100 percent thrust 
condition. The reflection factors are equal with 
magnitude Rj = R 0 = 0.35 and equal but varying 
phase angles. The variation of the phase angle, 

* = 0 through 2n radians, causes ilarge changes 
in the spectral shape. All spectra are double 
humped in the frequency range shown. It appears 
that more exact knowledge of the phase angle is 
required if one expects to accurately predict the 
effect of reflections on combustion generated noise 
spectra as obtained from combustor pressure 
measurements. 

Comparison of Combustor Data With Theory 

The narrowband fluctuating pressure spectra 
from Ref. 7 are shown in Fig. 6 along with the 
spectra predicted by the present theory. A range 
of operating conditions from 99.9 percent to . 

3.8 percent of design thrust are shown, (Fig. 6(a) 
through 6(d), respectively). As discussed in 


Ref. 8 the measured spectra above idle setting are 
bi-lobed with maximums occurring in frequency 
regions of 350 and 1100 to 1300 Hz. The theoreti- 
cal spectra were calculated for Rj = R 0 = 0.35 and 
a = n/2 and are shown in Fig. 6 by the solid lines. 
Also shown in Fig. 6 are the calculated spectra 
for reflection factor equal zero (dashed line). 
Reference 29 indicates that reflection factors on 
the order of 0.4 might be expected from turbines. 

The theory predicts two humps in the spectra, 
one centered in the frequency range between 200 
and. 700 Hz and the other centered in the frequency 
range 1000 to 1600 Hz. The predicted humps, above 
idle condition, appears to occur at a higher fre- 
quency than the humps in the data by 100 to 200 Hz. 
Comparisons of the theoretical hump center fre- 
quency with increasing engine speed above idle 
shows that it is predicted to increase with speed. 
The experimental data also show increasing hump 
center frequency with increasing engine speed. 

The shape of the theoretical spectra and data also 
show similar trends with frequency. It has been 
shown that reflection factor phase angles and, at 
the higher frequencies, the measurement location 
have significant effect on the spectral shape. 

Even though the real phase angles and measurement 
location are not known precisely it appears that 
the agreement between theory and experiment is 
adequate. Taking into account the stated uncer- 
tainties the conclusion is drawn that the deviation 
of the spectral data above idle conditions from the 
simple source noise given in Ref. 27 may be attrib- 
uted to reflections from the turbine and combustor 
inlet. 

At idle condition (Fig. 6(d)), the theoretical 
pressure spectrum with reflection factors of zero 
(see dashed line Fig. 6(d)) agree with the experi- 
mental data. The data does not show a distinct 
peak in the spectrum where it might be expected if 
reflections were significant. Reference 29 indi- 
cates that the turbine reflection should decrease 
with turbine pressure ratio. Hence, at the 
3.8 percent thrust level the turbine reflection 
factor should approach zero. Judging from the 
lack of humps in the spectra at the 3.8 percent 
thrust level it appears that little or no reflec- 
tion exists in the combustor at idle speed and 
indeed the turbine reflection factor is zero. As 
has been shown in Ref. 27 the infinite tube theory 
predicts the narrowband spectrum in the CF6-50 
combustor very well at the idle condition. 

Concluding Remarks 

The present analysis neglects the effects on 
the acoustic pressure measured between the re- 
flecting planes of both area variation between 
reflecting planes and acoustic wave number varia- 
tion within the combustion zone. The analysis in 
Refs. 25 and 26 includes these effects. However, 
the present analysis does include the effect of a 
distributed source. Certain conclusions maystill 
be drawn from the present analysis that appear to 
explain some of the variations in the spectra from 
the simple combustion noise source derived in 
Ref. 27. 

It has been shown that reflections from the 
turbine and combustor inlets are the probable cause 
for the deviation of the combustor fluctuating 
pressure spectra from the smooth curve given by 
the source spectrum. Since the strongest reflected 
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wave interference with the source can only exist 
between the reflecting planes one may conclude that 
the source wave itself will dominate and essen- 
tially be transmitted through the reflecting planes 
without spectral shape modification. Of course the 
magnitude of the source wave will be decreased but 
the frequency dependence should remain the same. 

For this reason one may expect that the acoustic 
pressure measured external to an engine in the far 
field will have the very low frequency dominated 
spectrum of the combustion noise source as reported 
in Ref. 27. Cross correlations between combustor 
and far-field microphones verify this conclusion.® 

Conclusions 

A theoretical investigation of the effect on 
combustion chamber fluctuating pressures of re- 
flecting planes at the combustor and turbine inlet 
locations has been performed. The pressure spec- 
trum has been calculated and compared to existing 
data from a CF6-50 turbofan engine combustor and 
found to be in substantial agreement. The follow- 
ing conclusions have been drawn from the theory 
and comparisons of the theory with data. 

1. At idle engine speed no reflection from 
the turbine or combustor inlet occur and the in- 
finite tube theory applies and yields excellent 
agreement with the data. 

2. Above engine idle conditions, reflections 
from the turbine and combustor inlets occur and 
reasonable agreement between theory and narrowband 
combustor pressure spectra has been found using a 
reflection factor on the order of 0.35 and a phase 
angle of »'/ 2 radians. 

3. Measured spectrum shape is independent of 
the measurement location at low frequencies, but 
not at the high frequencies. 


L q distance to upstream reflecting plane 

from point source, m 

Lj distance to downstream reflecting plane 

from point source, m 

SPL sound pressure level, dB (re 20 p Pa) 

P instantaneous or mean static pressure. Pa 

p time-fluctuating part of instantaneous 

pressure. Pa 

P ref reference pressure used in calculating 

sound pressure level, 20 g Pa 
P u Fourier-transformed fluctuating pressure, 

Pa sec 

R reflection factor 

T mean temperature, K 

T' temperature fluctuation 

V mean velocity, m/sec 

v velocity fluctuation 

W mean mass flow rate, kg/sec 

x axial distance measured from combustor 

inlet, positive in direction of flow, m 
x Q point source location, measured from 

upstream reflection plane, m 
« time-fluctuating component of density, 

kg/m 3 

n c combustion efficiency 

n p acoustic pressure efficiency 

* phase angle, deg 

x fuel mass decay constant, assigned a 


Symbols 

A 

C 


P 

C 1’C 2 



Appendix A 


2 

cross-sectional flow area, m 
sonic velocity, m/sec; constant of 
integration 

specific heat at constant pressure, 0/kg K 
constants used in complementary function 
for solving differential equations 
functional notation, Eq. (B13) 
frequency, Hz; function 
fuel/air ratio 

2 

acceleration of gravity, m/sec 
heating value of fuel, J/kg (cal/kg) 

acoustic wave number, U /C, 1/m 
distance between reflecting planes, m 
burning length (i.e., axial distance 
measured from flame to point in 
combustor where chemical reaction 
ends), m 


value of 2 h/L k 

D 3 

p mean density of fluid, kg/m 

a Fourier-transformed source terms 

u 

u angular frequency, 2nf, rad/sec 

Subscripts: 

b burning length 

L left running 

R right running 

T total conditions, turbine location 

u Fourier-transformed quantity; spectrum 

0 upstream axial location 

1,T axial location downstream of combustion 

zone 

3 third source term (Ref. 27) 

Appendix B 

Solution of the One-Dimensional, Fourier Trans- 
formed Acoustic Wave Equation with Combustion Noise 
Source Term and two Reflecting Planes ’ ~ 

The acoustic wave equation governing the 
acoustic pressure generated and measured in a 
one-dimensional duct with reflecting planes 
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representing the combustor and turbine inlets is, 
from Ref. 27 : 


3 2 P 2 

— r + k p u = " 9 3 ,o 

ax 


(Bl) 


U « - X (x-L„) 


W u , o"c"p V 0 


9 3, u = -*r-y^- 


1 (B2 j 


p = C 

K tu,0 


r ik o x „ - i <V"*o ) T A 

1 [ e + R 0 e J, 0 < x < x 0 


< 8 . 
(B7 ) 


where the source term <?, is given in Ref. 27 as: 

j,<o 


In the region downstream of the point source the 
acoustic pressure again written in terms of the 
turbine reflection factor is: 
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The homogeneous equation obtained by setting 
to zero is 
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The solution of equation (B3) is 
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Wave Pressure Wave Pressure Wave Pressure 


A point or planar source is assumed to exist in 
the combustion zone. The regions between the point 
source and the combustor inlet and the point source 
and the turbine inlet have different temperatures 
and hence different sonic velocities and fluid 
properties. At the combustor inlet reflecting 
plane (x = 0) the reflection factor, a complex 
number., defined as the complex ratio of the reflec- 
ted to the incident acoustic pressures, is given 
by: 




(B5) 


and at the turbine reflecting plane, (x = i). 




(B6) 


The magnitude and phase of both reflection factors 
are assumed to be constant and known quantities. 

The two regions are treated separately until 
later in the derivation. In the region upstream of 
the point source equation ( B4 ) may be written in 
terms of the combustor inlet reflection factor as: 


The method used in Ref. 28 for solving the harmoni- 
cally driven string with reflections is now applied 
to find the solution of the acoustic wave equation 
used herein. Equations (B7) and (B8) have two 
unknown constants, Cj and C?. Therefore two bits 
of information are required to solve for these con- 
stants. The first requirement is that the pres- 
sures across the point source are continuous, i.e., 
pressures at the source are equal for both regions. 
Equating equation (B7) to equation (B8) at the 
source location, x = x 0 yields 
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From Ref. 28 the difference in the slopes at the 
point source is minus one. Taking the derivative 
with respect to x of Eqs. (B7) and (B8), equating 
the differences of the slopes to -1 yields the 
second expression in terms of Cj and C? as: 
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Solving Eqs. (B9) and (BIO) for Cj and C 2 and 
substituting and C 2 into Eqs. (B7) and (B8) 

respectively yields the expression for the acoustic 
pressure that would be measured in the duct if the 
duct were forced by the unit impulse function. The 
resulting acoustic pressures for the two regions 
are: 
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where the denominator is a function of the point 
source location and is given by: 
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Having the solutions for the acoustic pressure, in 
their respective regions as a response to the unit 
impulse function the acoustic pressure in the duct 
as a response to the distributed source in the com- 
bustion region can be found by integrating the 
product of the source 93 and the pressure 
response to the unit impulse function over the com- 
bustion region as discussed in Ref. 28. In equa- 
tion form the pressure in the upstream region is: 
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In the downstream region: 
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where F 0 (x 0 ) is given by equation (B13) and the 
source term from Ref. 27 is given by 
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The integration is performed by converting from 
the complex polar form to the complex rectangular 
form, integrating the real and imaginary parts 
separately, and then converting the results of the 
integration back to polar form. Numerical integra- 
tion was used in the present work. The pressure 
calculated using Eq. (B15) and expression (B16) is 
the theoretical acoustic pressure in the duct for 
the region downstream of the point source. It was 
used in making. the theoretical acoustic pressure 
predictions reported herein since the data used 
for comparison was taken from a sensor located 
downstream of the expected flame front location. 
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Figure 3. - The effect of measurement location on acoustic 
pressure spectrum calculated from theory: 99.9-percent 
thrust, CF6-50 turbofan operating condition, equal phase 
angles (0 = tt/ 2) and reflection factors (Rif 0. 35), L g/£ = 
niB i_i» = n 89 i.i* = i nn » = n w m n=n nim 
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(b) Reflection factor of 1.0. 

Figure 4. - The theoretical effect of reflection factor magnitude 
on spectral shape for CF6-50 combustor operating at 99. 9- 
percent thrust condition; equal reflection factors and phase 
angles (O = irl 2 rad. ), L Q li - 0. 18, L T /£ = 0. 82, L b /i = 1. 00, 
£ = 0. 348 m, x/i = 0.2, rip' 0.883. 
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Figure 5. - The theoretical effect of phase angle on spectral 
shape for the CF6-50 combustor operating at 99. 9-percent 
thrust level; equal reflection factors (R = 0. 35) and phase 
angles; ■ 0. 18, L T ll ■ 0. 82. U/£ = 1. 00, x/i = 0. 20, 
i » 0. 348 m, rip = 0. 08^3. 
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